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ABSTRACT 



The effective equations of motion for low-frequency mean gauge fields in the Abelian Higgs 
model are investigated in the presence of a scalar condensate, near the high temperature equi- 
librium. We determine the current induced by an inhomogeneous background gauge field in the 
linear response approximation up to 0(e 4 ), assuming adiabatic variation of the scalar fields. 
The physical degrees of freedom are found and a physical gauge choice for the numerical study 
of the combined Higgs+gauge evolution is proposed. 
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1 Introduction 



The real-time transformation of light gauge fields into massive intermediate vector bosons is a 
subject of increasing cosmological interest. 

Recently, Garcia-Bellido et al. [l[ (see also [§]) assumed that the coupled inflaton+Higgs 
system was far out of equilibrium when the energy density of the expanding Universe has passed 
the point corresponding to the thermal Higgs transition. With the supplementary condition 
that the reheating temperature after the exit from the inflationary period did not exceed the 
electroweak critical temperature, they demonstrated in a (1+1) -dimensional toy model, that 
large enough matter-antimatter asymmetry could have survived till today. The existence of 
non-equilibrium Higgs transitions in 3+1 dimensions has been demonstrated by Rajantie et al. 
[Q. In these investigations the classical mean field equations were used with renormalized and 
temperature-dependent couplings. 

A more traditional field of the investigation is the evolution of the baryon asymmetry 
through the electroweak phase transition, accompanying the onset of the Higgs effect. It has 
been thoroughly discussed, under the assumption that the system is in thermal equilibrium, 
with specific emphasis on the high-temperature sphaleron rate 0, |5|, ||. Near equilibrium, the 
gauge-invariant HTL action dominates the influence of the thermalized quantum fluctuations 
on the motion of the mean fields with k ~ O(gT) ]7], |8|, [|. This term was taken into account 
in real-time simulations of the sphaleron rate || (see also ||10||). 



The framework for the theoretical study of the baryon asymmetry has been somewhat 
modified with the realization that, within the Standard Model, the Higgs effect sets in via 
smooth phase transformation, characterized by an analytic variation of the order parameter 



IJ], [12], [L3[ [14]. Under such circumstances the expectation value of the Higgs field is different 
from zero also above the electroweak energy scale. In this context it is remarkable that the 
leading HTL correction is insensitive to the presence of the scalar condensate IfO) , [3] . 

In this note we go one step beyond the derivation of the HTL correction. Our aim is to 
find the leading effect due to the presence of an arbitrary constant background field: $. In 
equilibrium, such a background is sufficient for the calculation of the effective potential, but not 
the full effective action. Similarly here, we restrict somewhat the generality of our correction 
terms to the equations of motion, namely the corrections to the $ 2 -A 2 vertex will be determined 
neglecting the non-local effects in the scalar field. 

The high-temperature limit of the resulting expression of the induced current shows that the 
next-to-leading "mass and vertex" corrections can be uniquely decomposed into gauge-fixing 
invariant and gauge-fixing dependent modes. In this way we can propose an expression of wider 
applicability for the Landau damping effect in the presence of a scalar condensate. 



Baacke and his collaborators |16l [L7L 18 1 have applied, in a series of papers, a complemen- 



tary approach to the Higgs effective action without assuming the presence of a general gauge 
background. 

A non-trivial mean gauge field configuration reacts back on the scalar condensate. Our final 
goal is to derive a coupled set of equations for the scalar condensate and the mean gauge field, 
which is appropriate for studying the real-time onset of the Higgs effect and the variation of the 
sphaleron generation and decay rate. We propose a physical gauge for the numerical solution 
of these equations where the gauge-fixing dependent mode does not propagate. 
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The generalization of the derivation to non-Abelian systems presents only technical compli- 
cations. The Abelian calculation has been performed both in the standard real-time generator 
functional formulation and with an iterated solution of the linearized Heisenberg equations 
of the quantum fluctuations [[ 



21, 22 



We shall use below the latter approach, where the 
adequate two-point functions enter more naturally. The general linear response theory will be 
described in a separate publication |23| . 



2 Mean Field Equations in the Abelian Higgs model 

The Lagrangian of the model in the 0(2) notation is the following: 

£ = —F^Ff + \(d^f - \m 2 & + ei M > + ^i 2 $ 2 - A(l> 2 ) 2 , (1) 

where <S = ($i, $ 2 ) and J M = $ 2 <9 M <&i - <f>i<9 M $ 2 . We split the fields into a mean field (A^ $) 
and a fluctuation contribution 

A M = + a M , <f> = $ + <p, (a„) = ((p) = (2) 

at any time (averaging is understood with respect to the initial density matrix). We assume 
that the scalar background is constant and that it points to the $i direction, its value being $. 
We use in the sequel the notations mw = e$, m 2 H = m 2 + Ai> 2 /2, m 2 G = m 2 + A$ 2 /6, although 
these are not the vacuum values. 

We fix a 't Hooft gauge by changing the Lagrangian as 

C^C- —{d^ + Zm W (p 2 ) 2 - c(d 2 + im 2 w + e£m wVl )c, (3) 

where c is the ghost field. In this way the field (p 2 receives an additional mass contribution 
im 2 w . 

The operatorial equation of motion (EOM) separately induces EOMs for the mean fields 
and the fluctuations. The average of the operatorial EOM for the gauge field reads as 

A » + fi«(A) = 0, (4) 

with the induced current 

f; d (A) = e M > + 2e 2 $ ( W > + e 2 A /I (<p 2 ) + e 2 (a^ 2 ) , (5) 

where j M = ^p 2 d^\ — ( Pi9/ J ,<p 2 - We will perform the calculations at the one-loop level, when 
the last term does not contribute. The subtraction of (|]) from the full equation yields the 
EOM for the fluctuations. At one loop it is sufficient to consider only the equations linearized 
in the fluctuations. On the other hand, since we want to calculate in the linear response 



(d 2 + m 2 w )g. 
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approximation, we can neglect all terms non-linear in A. These assumptions make the equations 
very simple: 



1 



{d 2 + m 2 w )g^+\l--jd l ,d v 



a" - 2e 2 <S>A^i = 0, 



(<9 2 + m|>! + 2eA»d fM ip 2 + e{dA)^ 2 - 2e 2 $A l a tl = 0, 
(8 2 + m 2 G + £rn 2 w )y 2 - 2eA»d fM ip 1 - e{dA) Vl = 0. 



(6) 



The ghost fields follow a free EOM, so that they do not influence the present calculation. These 
equations are solved to linear order in the A background: 



aJx) = a<£\x)-2#* I *zG%{x - z)A"(z)tf> \z) 



R 



^i(ar) = + e ld 4 zG K (x - z) 2A"{z)d^ 2 , {z) + (dA)(z)^ 2 > (z) - 2e^A"(z)a^(z) 



(0), 



,(0), 



<p 2 {x) =$\x)-e l<?zG?{x-z) 2A»{z)d^r(z) + {dA){zW(z) 



(o), 



(7) 



where the superscript zero denotes the solutions of the free EOM, and G R, s are the free retarded 
Green functions^. 

For the induced current (^) we need the expectation values of certain local products of the 
fluctuating fields. For example, it directly follows from the above equations (with (AB) = 
that 



^ 2 {x)d^ x {x)) = eld 4 z\ d£G*(x-z) 



2A v {z) fa{x)&l<p*(z)) Q + {dA){z) (Mx)Mz)) 



- G?(x-z) 



2A\z) {d lt <p 1 (x)d u ( P i(z)) + (dA)(z) (d^iWpiiz)), 



(8) 



We define the local products in a symmetric way (i.e. \ ((p 2 (x)d^(fii(x) + d fl <fi(x)ip 2 (x))) and 



introduce 



A(x -z) = - ((f(x)(f{z) + ip(z)ip(x)) . 



(9) 



Then the above expression can be written in Fourier space as 

d A p 



fadpipi) (Q) = -eA"(Q) 



P,(2p-Q) u G H (p)A 2 (Q-p) + G R (Q-p)A 1 (p) . (10) 



{2-kY 



The evaluation of other expectation values goes along the same line, finally giving 

d A p 



e{j,) (Q) = -e 2 A»(Q) 



(2tt)< 



(2p - Q)»(2p - Q) v G R (p)A 2 (Q -p) + G R (Q - p)A 1 (p) 



2e 2 $ ( W > (Q) = -4e 2 m 2 w A»(Q) 

d 4 p 



d 4 p 



G R {j>)A x {Q -p) + G R (Q - p)A^{p) 



e 2 A, ( v 2 ) = e 2 A,(Q) J ^JL [Al (p) + A 2 (p)\ . 



'IF 



4 Here the definition KG R = —6 is used, where K denotes the free kernel. 
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Assuming that the free fluctuations are in thermal equilibrium, the propagators can be related 



to the corresponding spectral functions ]24], which are the discontinuities of the free kernels of 
eq. (H). Introducing 

1 



A m2 (p) = 2n(^ + n(\p \)^6(p 2 -m 2 ), and G^{p) 



p2 _ m 2 _|_ i £ p 



where n is the Bose-Einstein distribution, we have 



Ai = A 



A 2 - A m 2 aHm 2 w , 



A 



-9p,vA m 2 w + (A r 



m 



,2 



(12) 



(13) 



and analogously for the corresponding G R, s. 



3 High-temperature expansion 

The leading HTL term of j™ d in the high-temperature expansion comes from e by neglecting 
all the masses [15]. Our aim is to calculate the first subleading term proportional to $ 2 instead 
of T 2 in the high-temperature expansion. 

We emphasize that j (Q) = j md {—Q)*, therefore all partial contributions that are odd 
under hermitian Q-reflection can be freely omitted. 

A useful formula, which simplifies the calculations, is the following [25|| : 

I F (Q,m\ M 2 ) = J -0- F(Q,p) (gUp)&aMQ-p) + GUQ - P)A„ 



(2tt) 4 2 ^' 2pQ + m 2 -M 2 

The tensorial structures appearing in ( |TTD when cast into the form (|T4] ) imply the appearance 
in F of the p-dependent terms Pfj,Pu,P^Qu + Pi/Q/x anc ^ a ^ so OI " terms independent of p. 

3.1 Even terms 

We start the discussion with the terms even under p-reflection and introduce the notation 
f(p) = (F(Q,Q/2 — p) + F(Q, Q/2 + p))/2. Here the m 2 — M 2 term in the denominator can 
be neglected since in the mass expansion 



(2jt)"^" ' (2pQ)= 

because of the odd p —>■ —p behaviour of the integrand. Then we only need to expand the 
difference of A's of the numerator. The numerator can be expanded with respect to Q, when 
low-momentum mean fields are considered. The leading term of the gradient expansion gives 
zero, again because of the p-odd integrand. 
The first non-zero contribution is therefore 

/,(„», m 2 ) = 4/1^ M|_ (An>1 w + Ami(p)) . (16) 
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Now, we treat separately the two actual cases: / = p^p v and / = constant. The / = p^p v 
case contributes to the full expression of j™ d {Q) 

2e^( W /(g; H£(A, + A 2 + " A^). (17) 

We introduce the field-strength tensor by A U Q P = —iF vp + A P Q V . The local term proportional 
to A reads as 

d A p d 



e AP{Q) I (#^ (Al + As + " A ^ } - (18) 

After partial integration the first two terms cancel with e 2 A fl (ip 2 ) in the induced current. What 
remains is a contribution vanishing at zero m^. In the mass expansion they give 



■local,l 
J fj, 



^-O^A,(Q)J^-^. (19) 



(2vr) 4 dp 2 ' 

For each A m 2 the field strength contribution of fll7D is rewritten with the help of the relation 

„<9A dA dn(\pn\) n <^ , s 

s V = - 2ps s? +So ^ 2rf(p m) (20) 

in the form 

- - 2 /l^ + fw^gW - »■>) ■ PD 

The first term drops out because of the antisymmetry of F. In the second term we perform the 
mass expansion. After adding the different contributions we find 

d 4 p PyPv dn(\p Q \) 
(2tt) 4 2pQ dp 

The first term of this expression is the usual HTL contribution, the second is a mass correction 
to it. In the mass correction, originally, there was m 2 + m 2 + — ^m^y, but because of 
m i = m G + i m w the £ dependence drops out. 

The / = constant contribution to the induced current appears as 

d A p 1 dAp 
(2tt) 4 2pQ dpp 

Since it is already proportional to we have dropped the mass dependence of the integral. 
Using once more the relation ( PUD we write for it 



(i) = 4e 2 MF I/0 ' 



n^A". (22) 



The first term is again local: 



^ = 4e=m^(Q)/ I g I ^. (25) 



The second term reads as 
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3.2 Odd contributions 

Odd contributions come exclusively from (a^i). Using ( |14|) and denoting the current by jjp 
we find 

Finally, with the help of ( p0|) , performing the mass and external momentum expansion according 
to the method followed in the previous subsection we arrive at 

3.3 Linear response and physical modes 

The full induced current is the sum of the different parts coming from (|i~9[) , (^2|) , (|25|) , (p6|) and 



(i 



W = J ^ +J (l) +J (2) +J (3) ; (2g) 



^local ^local,l i Aocal,2 

Ifj, 



^ ^ 3 J fj, ^ «7 ^, 

Here ffl and j^ 2 -* are ^-independent, while j l ° cal and depend on the gauge fixing. For 
the physical characterization of the system (for instance, damping rates) we have to find the 
independently evolving modes. In order to do this we decompose the polarization tensor in the 
tensor basis, appropriate for finite-temperature studies 



Q ^5 Q 

(30) 

where Q = (q , q), q 2 = q 2 and ti^ = g^ - Q^qo/Q 2 was used. 

Since of fl22|) is transverse (Q^U^ = 0), it is the combination of P T and P L . Introducing 
n L = IToo and n T = l/2P^ v W v we find 

U = -^-U L P L + U T P T . (31) 

The Z^y term in ( PH[ ) can be written as a combination of P G and S. Introducing Z G = Zjj and 
Z s = (q 2 Z G - Q 2 Z 00 )/(qq ) we find 

Z = Z G P G + Z S S. (32) 

Using the completeness relation P L + P T + P G = g the Fourier transform of the EOM (^) has 
the following form: 



(Q 2 - R - U T )P T + (Q 2 - R + |V) P L +(^Q 2 -R-Z G )P G -Z S S 
with R defined from m^A^ + j l ° cal + j^ 2 ) = iL4 M . 



A = 0, (33) 



Since S is not a projector, but mixes the subspaces belonging to P L and P G , in this two- 
dimensional subspace the inverse propagator matrix still has to be diagonalized: 



-Q 2 — R — Z, 



G 



(34) 



The eigenvalues should be found with one-loop accuracy, which means that terms proportional 
to the square of one-loop corrections (i.e. U 2 L , Zq and Zg) should be neglected. This, however, 



implies that the eigenvalues are the diagonal entries, 
the (P G , P L ) plane, we can therefore write 



With appropriately rotated projectors in 



(Q 2 -R- U T )P 



T 



Q 2 




J G 



P 



G 



A = 0, 



(35) 



As expected, the transverse modes plus the Pl mode, which might be called longitudinal, can 
be made independent of the gauge-fixing parameter (see the discussion on the renormalization 
of R below). 



3.4 The infrared separation scale 



The p-integrals in all terms of j™ are factorized into a radial and an angular integral. It is well 
known that the HTL term describes the dynamical screening of the electric fluctuations below 
the Debye scale, an infrared separation scale ( "IR cut-off" ) has to be therefore introduced into 
the radial integration at p = M = Cm x eT. The effective equations one arrives at in this way 
are to be used for the modes p < M. 

Applying this cut-off to the integrals appearing in the expressions of R(Q) and H^ u , one 
finds 

R = 



n, 



IT 




+ m G + m 



H 



go ln go + g 
g go - g' 



(36) 



with k = 27rexp(— je) and m 2 D = e 2 T 2 /3 — e 2 MT/ir 2 . The logarithmical UV divergence in R 
can be absorbed into the e 2 renormalization. With appropriate renormalization scale fi = kT 
the £ dependence can be made to vanish. 

The intermediate IR scale M contributes a term to the Debye mass, which depends linearly 
on M. It will be cancelled by the linear M-dependence of the self-energies, which shows up in 
the (one-loop) solution of the classical effective EOM Qf) gjj, 0, ]30|. 

The mass corrections in (|35|), proportional to & 2 A U , should be interpreted as coming from 
an induced (non-local) $ 2 A 2 vertex correction. In the subsequent classical time evolution it 
contributes to the self-energies of both fields. When (perturbatively) combined with the ~ MT 
terms coming from the tadpoles, it results in a finite contribution of order T 2 , independent of 
the exact choice of the coefficient Cm in the expression of M. Therefore, these terms can by 
no means be neglected. 
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4 Conclusions 



In this letter we have investigated the induced current in the effective EOM of the gauge 
field in the Abelian Higgs model. Beyond the well-known HTL term, we now determined also 
the subleading ~ ^A^ corrections in adiabatic approximation ($ = constant). The most 
important property of the corrections is that, just as the leading HTL term, in the high- 
temperature expansion for the physical degrees of freedom it is independent of the gauge-fixing 
procedure (after appropriate renormalization). 

For consistency we have applied an IR cut-off M = Cm x gT to the fluctuations, and the 
effective EOM is valid below this scale. The subsequent 3D time evolution will be insensitive to 
the accurate choice of M. Partly it is cancelled by the 3D (Rayleigh- Jeans-type) divergences, 
partly it yields also non-zero, M-independent contributions, when the gauge-scalar vertex 
corrections are combined with the 3D would-be divergences. 

We shortly discuss here the numerical implementation of the effective gauge field equation. 
The most convenient gauge fixing seems to be the Landau-gauge, where the ^-dependent mode 
does not propagate. Then can be left out of the discussion and the effective equations of 
motion are to be solved under the constraint <9 M A M = 0. It can be implemented, for example, 
by solving the equations for Ai and computing Aq from the constraint. 

The non-local induced currents are written in local form with the help of auxiliary fields. 
With the well-known form of the leading HTL current, we can write in this form also the 
corrections due to the non-zero scalar background: 

, (1) 2 r dtt v ( e 2 T m 2 w + m 2 G + m 2 H \ w 
3\ = m D J [v iVj - — Sij j Wj (x, v), 

.(2)/ \ e 2 m 2 v T rdfl v (2), , 

3i {x) = -^rJ (x,v) ' (37) 

where the auxiliary fields satisfy 

(d - vd)W^ = F 0l , (d - vd)wl 2) = d A. (38) 

In the derivation of the expression of we have performed a partial p-integration in the 
part of (|22|) proportional to S'(p 2 ), and used the fact that, to the accuracy of our calculation, 
Q u F utJ- (Q) = can be exploited in the expression of the induced current. 

The equations for the scalar field receive mainly local corrections (renormalization and In- 
dependence of the couplings in the classical equations). However, for consistency also the 
^-derivative of the $ 2 -v4 2 vertex correction established in the present calculation should be 
introduced. 

The logics of the derivation followed in the Abelian Higgs model seem to be robust enough 
for us to attempt its generalization to the non-Abelian case, which will be the subject of a 
future study. A numerical study of the corrected EOMs can lead us to a deeper understanding 
of the reliability of our present views on the high-temperature Higgs models, especially where 
the sphaleron rate is concerned. We will also learn in more detail of the non-equilibrium aspects 
of the cosmological onset of the Higgs regime. 
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